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Abstract 

Inclusive K° and A photoproduction has been investigated at HERA with the 
HI detector at an average photon-proton center of mass energy of 200 GeV in 
the transverse momentum range 0.5 < pt < 5 GeV. The production rates as a 
function of pt and center of mass rapidity are compared to those obtained in 
deep inelastic scattering at (Q 2 ) = 23 GeV 2 . A similar comparison is made of 
the rapidity spectra of charged particles. The rate of strangeness photoproduc- 
tion is compared with pp measurements. The observations are also compared 
with next-to-leading order QCD calculations and the predictions of a Monte 
Carlo model. 
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1 Introduction 



The HERA collider allows the study of ep interactions over a wide range of squared four- 
momentum transfer, — Q 2 , from photoproduction (Q 2 ~ 0) to very high photon virtuality. 
This makes it possible to investigate the Q 2 dependence of particle production. 

Photoproduction interactions are in many respects very similar to hadron-hadron inter- 
actions JI| . As Q 2 increases, the ep interaction is considered to be deep inelastic scattering 
(DIS) in which the exchanged virtual photon interacts directly with a parton in the proton. 

The hadronic final state resulting from the fragmentation of partons is of a non-perturba- 
tive nature and is usually described using phenomenological models such as the string 
fragmentation scheme implemented in Jetset . In contrast to the study of unidentified 
charged particles || [|, strange particles allow tagging of a specific quark species. Strange 
quarks are mainly produced in the hadronization phase, which is dominated by particles 
with low transverse momentum, p t . However, they can also originate from hard partonic 
radiation and the direct interaction of a photon with the proton, e.g. boson-gluon fusion, 
thus contributing to the high-p t tail of strange particles. 

In this paper, the transverse momentum and rapidity spectra of K° and A particles 
in photoproduction at an average center of mass energy (W) = 200 GeV are presented. 
These measurements are made in the photon fragmentation region, 1.3 < y* < 2.8, and 
the transverse momentum range 0.5 < p t < 5 GeV and are compared with results obtained 
in DIS at (Q 2 ) = 23 GeV 2 and a similar average W. The comparison is extended to the 
rapidity spectrum of charged particles to see if the similarities observed for strange particles 
are of a more general nature. This allows the investigation of the Q 2 dependence of particle 
production rates in the photon fragmentation region. In the central rapidity region these 
rates are generally assumed to be independent of Q 2 [[J . The process dependence of particle 
production rates is studied by comparing the rates of strangeness production in jp and 
in pp. Finally, results from next-to-leading order (NLO) QCD calculations || and from 
Pythia |7], |^] are compared with the measurements. 

2 Models of Photoproduction 

The Pythia 5.7 Monte Carlo program, as used here, generates events according to the 
description given in ]7|, 0. The structure of the photon is parameterized according to |7|, [|. 
Furthermore, the Pythia option of generating multiple interactions within the same event 
is used 11. Fragmentation is performed using the Lund string fragmentation scheme, 
as implemented in Jetset 7.4. Matrix elements are calculated to leading order, and 
higher order terms are simulated by parton showers in the leading log approximation. The 
production rate of strange particles in Jetset is mainly controlled by the strangeness 
suppression factor, X s , that is, the probability of producing a strange quark pair relative 
to a light quark pair. By default X s is set to 0.3. However, recent deep inelastic muon 
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nucleoli scattering data from E665 | ID| , e + e data from DELPHI [[□]] and ep data from 
HI and ZEUS @ are better described by Jetset if X s ~ 0.2 Q. 

The authors of calculate the inclusive single-particle photoproduction cross sections 
as follows. The hard subprocess is calculated to next-to-leading order in QCD, as is the 
evolution of the fragmentation and parton distribution functions. The fragmentation func- 
tions are derived from e + e~ collider data. The parton distributions for the photon and 
the proton constitute an important input to the cross-section calculations. The numerical 
calculations presented here are based on the GRV parameterizations given in fl~3||. 



3 Experimental Procedure 
3.1 Detector 

The results presented in this paper are based on data taken in 1994 with the HI detector. 
During this running period 820 GeV protons and 27.5 GeV positrons were brought into 
collision in HERA. 



A full description of the HI detector can be found in |14|. In this paper only the 
components of relevance to the measurements presented here are mentioned. These are 
the positron tagger, the central tracker, the backward proportional chamber (BPC), the 
backward electromagnetic calorimeter (BEMC) and the liquid argon (LAr) calorimeter. 
The direction of the z-axis is chosen to be along the proton beam direction. The polar angle 
9 is defined with respect to the z-axis and the pseudorapidity is given by r\ = — ln(tan6'/2). 

The positron tagger, located 33 m from the interaction point in the positron direction, 
measures the energy of the scattered positron and, in conjunction with a photon detector, 
the luminosity by exploiting the Bethe-Heitler process [ 15 ]. 



The central track detector consists of an inner (20 < R < 45 cm) and an outer (53 < 
R < 85 cm) concentric central jet chamber (CJC), multi-wire proportional chambers for 
triggering purposes and two additional drift chambers which measure accurately the z- 
coordinate. The pseudorapidity range covered by the central track detector is \r]\ < 1.5. 

The LAr calorimeter, which surrounds the tracking system, consists of an electromag- 
netic and a hadronic section. Here it is only used to identify and reject events containing 
a large rapidity gap in the forward direction, in the range 2.03 < r] < 3.26. 

In the backward direction the BPC is used to determine the polar angle and the BEMC 
the energy of the scattered positron in DIS events. The pseudorapidity range covered by 
these detectors is —3.35 < r\ < —1.51. 

1 The following Jetset hadronization parameters have been changed with respect to the 
default settings (default/DELPHI/E665): Parj(2)=A s =(0.3/0.23/0.2), Parj(ll)=(0.5/0.365/0.5), 
Parj(12)=(0.6/0.41/0.6). Parj(2) is the strangeness suppression factor, Parj(ll) is the probability that 
a meson containing u or d quarks has spin 1, Parj(12) the probability that a meson containing s quarks 
has spin 1. 
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The calorimeters are surrounded by a superconducting coil providing a uniform magnetic 
field of 1.15 T in the region occupied by the central tracker. 



3.2 Event Selection 

Photoproduction events are selected by requiring the presence of a scattered positron in 
the small angle positron tagger. For these events Q 2 < 10~ 2 GeV 2 . The inelasticity of the 
scattering is given by y = 1 — |^ ■ sin 2 9/2. Here, E e is the energy of the incident positron, 
E' e and 6 the energy and the polar angle of the scattered positron. In photoproduction, 

E' 

to a very good approximation, the inelasticity is given by y = 1 — -et. Photoproduction 
events are required to lie in the range 0.3 < y < 0.7, in order to confine the scattered 
positron to a region in which the acceptance of the positron tagger is well understood. The 
trigger for these events requires a signal in the positron tagger and the presence of at least 
one negative track in the CJC More information on this track-based trigger can be found 
in Q. 

The trigger and the event selection for the DIS sample are as described in except 



that the y range is restricted to ensure that the mean hadronic center of mass energies 
of the DIS and photoproduction samples are similar. The scattered positron is detected 
in the BEMC and is required to have an energy greater than 12 GeV and a polar angle 
between 150° and 173°. The kinematic variables are determined using the measured energy 
and angle of the scattered positron. The DIS sample is restricted to and corrected in the 
kinematic region given by 10 < Q 2 < 70 GeV 2 , 10~ 4 < x < 10~ 2 and 0.3 < y < 0.6, where 
x is the Bjorken scaling variable. 

For both the jp and DIS samples, the corresponding integrated luminosity is 1.3 pb _1 . 

In addition to the measurements of strangeness production in jp and DIS, results on 
charged-particle production in these two processes are presented. For this purpose, a 
charged-particle photoproduction sample is selected consisting of events triggered by the 
coincident detection of a positron in the small angle positron tagger and at least one track 
pointing to the interaction vertex with p t > 200 MeV measured in the cylindrical multi- 
wire proportional chambers. The charged-particle DIS sample consists of events which are 
triggered by a cluster in the BEMC with an energy deposit of more than 4 GeV. For both 
the DIS and 7p charged-particle samples, the inelasticity y of the interaction is restricted 
to 0.3 < y < 0.5, resulting in an average W of 187 GeV. The DIS sample is restricted to 
the range 8 < Q 2 < 30 GeV 2 giving an average (Q 2 ) = 15 GeV 2 . The charged-particle 
data samples studied correspond to an integrated luminosity of 1.5 pb _1 . 

Common to all four event samples is the requirement of a reconstructed primary vertex 
within 30 cm of the nominal vertex position. 

The background from photoproduction in the two DIS data samples is estimated to be 
less than 3%. 

For a comparison of particle production in different processes it is important to take 
into account the different contributions of low-mass diffractive events. While in the case 
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of DIS this fraction is only 10%, it is around 30% in photoproduction ]!? , I^fl . These 
events are characterized by the absence of energy deposited in the forward direction. To 
allow a direct comparison between DIS and 7p these large rapidity gap events are removed 
from the samples by requiring the presence of at least 500 MeV of energy in the forward 
region of the liquid argon calorimeter, 2.03 < r] < 3.26. The remaining samples are termed 
non-diffractive in the following. 

3.3 Kg and A Identification in Photoproduction 

Kg mesons and A baryons are identified through their decay channels 

Kg -> 7r + 7r 

A (A) — >■ pn~(pit + ) 

where the pion, the proton and the antiproton tracks are reconstructed in the CJC The 
tracks of the decay products are required to satisfy \r]\ < 1.5 and p t > 180 MeV, where p t is 
defined with respect to the beam axis. This restricts the acceptance to a region where the 
reconstruction efficiency is high. In addition, it is required that the radial positions ri> eg and 
r en d of the first and last measured point on the tracks fulfill r^ eg < 30 cm and r en ^ > 37.5 cm 
and that their radial length be at least 10 cm. Candidates for neutral particles decaying 
into two charged particles, hereafter called V° candidates, are searched for by performing 
a constrained fit to each pair of oppositely charged tracks. The fit demands that the 
tracks meet at a common secondary vertex and that the two particle momenta of the 
decay products transverse to the path of flight of the V° candidate be opposite. The V° 
candidate invariant masses are calculated assuming that the decay products are two pions 
for Kg candidates and a proton and a pion for A candidates. 

To reduce the combinatorial background and ensure good reconstruction of the decay, 
the radial distance d T from the primary to the secondary vertex is restricted to 2 < d r < 
18 cm for the K s and 3 < d r < 18 cm for the A candidates. Note that A decays with 
d r < 3 cm have a lower reconstruction probability due to the large asymmetry in the p t of 
the decay particles. 

The contamination from A decays in the K s sample is reduced to a negligible level 
by demanding that the transverse momentum of the decay particles with respect to the 
momentum of the V° candidate be greater than 120 MeV. Likewise, the K s contamination 
in the A sample is removed by excluding those A candidates for which the ir + n~ invariant 
mass, m n7T , falls in the window 0.48 < m W7T < 0.52 GeV. 

The applied cuts ensure that kinematic regions in which the reconstruction efficiencies 
for K s and A start to degrade are avoided. The large number of Kg candidates makes 
it possible to exclude in addition regions in azimuth, tp, of the CJC which have known 
imperfections (and where systematic effects are larger) and to reject K s topologies which 
are dominated by combinatorial background [O. 
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In order to ensure an optimal acceptance of the CJC for the V°'s, their pseudorapidity 
has been restricted to \tj\ < 1.3 and their transverse momentum to 0.5 < p t < 5 GeV for 
K s and to 0.6 < p t < 5 GeV for A. 

The final K s and A signals, after applying all the aforementioned cuts, are shown in 
Fig. |l|. Approximately 7 700 K s candidates and 3 600 A candidates remain. In the fits the 
signals are described by a superposition of two gaussians on a linear background, accounting 
for varying invariant mass resolutions for different decay topologies. On average the RMS 
width of the K° s signal is 9.2 MeV, the RMS width of the A signal is 2.9 MeV. The fitted 
peak positions are in agreement with the PDG K° and A masses [2(J within the systematic 



errors resulting from performing the above fits using various functional forms and from 
remaining uncertainties in the calibration of the tracking system. 

3.4 Corrections and Systematic Errors 

In order to obtain the inclusive K° and A photoproduction cross sections it is necessary 
to correct the observed numbers of Kg mesons and A baryons for the branching ratios 
into the observed decay channels, and for the acceptances and efficiencies of the various 
detector components. 

The dependence of the acceptance of the positron tagger on the inelasticity, y, is known 
from a measurement of the tagger acceptance for the Bethe-Heitler process ep — *> epy and 
Monte Carlo studies of the luminosity system and the HERA beam optics. More details 
can be found in pi . 



The efficiency of the track trigger is studied by using a reference event sample obtained 
with an independent trigger which does not make any requirement on the hadronic final 
state. On average it is found to be 85% and only weakly dependent on the transverse 
momentum of the V°. The systematic uncertainty associated with the trigger efficiency is 
less than 6%. This has been cross-checked by using a sample of simulated and reconstructed 
events. 

A visual check of the efficiency for reconstructing tracks satisfying the requirements 
described in the previous section reveals that it is consistent with the estimate of 98% 



obtained in [2!|. The systematic uncertainty assigned to this correction factor is 2%. 

The main correction results from the limited geometric acceptance of the CJC for V° 
decays and the cuts applied to select V° candidates explained in the previous section. 
The acceptance can be parameterized as a function of the pt and rj of the V° and ranges 
between 3% and 17% for Kg and between 4% and 25% for A decays. The cuts which 
have a large effect on the acceptance are the lower cuts on d r and p t of the decay tracks 
and, in the case of the Kg, the ip cut. The bin-dependent systematic uncertainty on the 
acceptance calculation is estimated by varying all cuts independently and re-evaluating the 
cross section. In a given bin the maximum variation is taken as a measure of the systematic 
error resulting from the acceptance correction. The bin-dependent systematic errors are 
listed in Tab. 1 and 2 and range between 2% and 10%. 
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Within the region defined by the above cuts the efficiency for reconstructing l^ 0, s is 
found to be 100% using a full detector simulation of photoproduction events. Since the K° 
analysis is restricted to the fully efficient regions of the CJC, which are well simulated, no 
additional correction is applied. For the reconstruction of A decays a detection efficiency 
of 97% is found. This leads to a systematic error of tfflo for K s and 1|% for A. 

The background resulting from non-ep interactions (proton wall and proton gas inter- 
actions) has been estimated from the proton pilot-bunch events, that is, events where the 
proton bunch does not have a colliding positron bunch partner. This background has been 
found to be negligible, but a systematic uncertainty of 1°% is assigned to its contribution 
to the cross section. 

The systematic uncertainty of the luminosity determination is 1.5% [p3fl . 

The overall bin-independent systematic error is calculated by adding the various contri- 
butions in quadrature and results in errors of ±8% for the K° cross section and lf% for 
the A cross section. 

For the charged-particle spectra, the dominant contributions to the systematic error, 
both in photoproduction and DIS, originate from the extrapolation to pt = (5%). The 
extrapolation was performed using a sample of fully simulated Monte Carlo events. The 
resulting error was estimated by varying the minimum p t of the particles both in Monte 
Carlo and data. In the case of the total charged-particle photoproduction sample an 
additional uncertainty in the losses of events with low charged-particle multiplicity (6%) 
is taken into account. 

4 Experimental Results 

Comparisons between the photoproduction data and the above-mentioned models and 
calculations are performed in the laboratory frame. Comparisons of strangeness production 
in photoproduction with that in DIS and pp interactions are made in the photon-proton and 
the proton-antiproton center of mass systems, using the non-diffractive data sample. The 
transverse momentum, p t , and the rapidity, y*, are defined with respect to the direction 
of the exchanged boson for DIS and with respect to the direction of the positron beam for 
7p. In the case of pp, pt and \y*\ are defined with respect to the beam axis. 

In Tab. p] and |2] the total 7p cross sections for K°, A and A production are presented 
as a function of p t and rapidity y*. The y* range has been restricted to 1.3 < y* < 2.8 
to make the measurement insensitive to the 9 acceptance of the CJC. The p t spectra fall 
steeply whereas there is almost no dependence on y* in the measured rapidity interval. 

Here and in what follows, all K s results are multiplied by 2 and therefore correspond to 
K° and K° production; no distinction is made between K° and K° and the combination 
of both is denoted by K°. Furthermore, it can be seen that the A cross section is, within 
errors, equal to the A cross section. Therefore they are combined and subsequent references 
to A always mean A + A. 
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For photoproduction the ep — > V°X cross section a ep is related to the jp — > V°X cross 
section a w (W) by the Weizsacker- Williams formula [24 



where F 



'ep - ^p(W) ■ F 

0.0136 is the flux factor integrated over 0.3 < y < 0.7 and Q 2 < 10~ 2 GeV 2 . 



a) K°: 





1 do 








A?/* dpf 


stat. 


syst. 


p t range 








[GeV] 


[nb/GeV 2 ] 




0.5 0.6 


25500 


1400 


500 


0.6 0.7 


16700 


900 


300 


0.7 0.8 


11300 


600 


200 


0.8 0.9 


7800 


400 


300 


0.9 1.0 


4900 


300 


200 


1.0 1.2 


3200 


100 


200 


1.2 1.4 


1600 


100 


50 


1.4 1.8 


570 


30 


30 


1.8 2.2 


190 


20 


10 


2.2 3.0 


38 


5 


3 


3.0 5.0 


3.6 


0.9 


0.4 



b) A: 


0.6 


0.9 


1900 


100 


100 




0.9 


1.2 


630 


30 


40 




1.2 


1.7 


190 


10 


20 




1.7 


3.0 


17 


2 


1 














c) A: 


0.6 


0.9 


1700 


100 


100 




0.9 


1.2 


600 


30 


30 




1.2 


1.7 


200 


10 


20 




1.7 


3.0 


15 


2 


1 



Table 1: Inclusive K° (a) and A/A (b,c) total jp production cross sections 1/Ay* ■ do/dp 2 
for 1.3 < y* < 2.8. The systematic errors given in the table do not include an overall 
systematic uncertainty of ±8% for K° and 1|% for A as explained in section 3.4. 



4.1 Comparison with Models and Calculations 

In this section the predictions of Pythia in combination with Jetset, and the results of 
next-to-leading order calculations, are confronted with the measurements reported in this 
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a) K°: 







da 

dy* 


stat. 


syst. 


v* 


range 


[nb] 




1.30 


1.55 


13500 


600 


700 


1.55 


1.80 


12700 


600 


400 


1.80 


2.05 


13200 


600 


300 


2.05 


2.30 


11800 


600 


500 


2.30 


2.55 


12600 


600 


400 


2.55 


2.80 


12200 


600 


400 



b) A: 



1.300 


1.675 


1900 


100 


100 


1.675 


2.050 


1600 


100 


100 


2.050 


2.425 


1600 


100 


300 


2.425 


2.800 


1500 


100 


100 



c) A: 



1.300 


1.675 


1900 


100 


100 


1.675 


2.050 


1600 


100 


100 


2.050 


2.425 


1300 


100 


200 


2.425 


2.800 


1400 


100 


100 



Table 2: Inclusive K° (a) and A/A (b,c) total •yp production cross sections da/dy* for 
0.5 < p t < 5 GeV in the case of K 0, s and 0.6 < p t < 5 GeV in the case of A's. The 
systematic errors given in the table do not include an overall systematic uncertainty of 
±8% for K° and 1|% for A as explained in section 3.4. 



paper. 

Figure shows a comparison of the measured K and A cross sections in p t and r] with the 
predictions of Pythia using standard parameter settings (A s = 0.3) as well as parameter 
sets obtained from fits to DELPHI and E665 data, as described in section 2. In the case 
of K° mesons (Fig. |2| a,b) the measured cross sections are in reasonable agreement with 
the Monte Carlo predictions using the DELPHI or E665 settings whereas the predictions 
with the standard settings seem to overestimate the cross section. However, the measured 
A (Fig. Q c,d) cross section lies significantly above the Monte Carlo prediction for all three 
scenarios. The UA5 collaboration has reported a similar discrepancy between A production 
in pp interactions at 200 GeV and Pythia predictions P5j . No attempt is made here to 
tune the Jetset parameters controlling diquark production. 

Recent next-to-leading order QCD calculations of inclusive particle spectra in ep scat- 
tering use fragmentation functions which have been fitted to e + e~ data 1]]. These NLO 
calculations involve four scales: the factorization scales connected with the parton densities 
of the photon and the proton, the factorization scale of the fragmentation functions and 
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the renormalization scale of the QCD coupling constant. The calculations for K° cross 
sections assume that all four scales are equal; they are taken to be equal to pt/2, p t or 
2p t (where p t is the transverse momentum of the K°). This arbitrariness in the choice 
of the scales introduces a change in the predicted cross sections and is a measure of the 
contribution of higher-order terms. To reduce the scale influence, the calculation for the 
photoproduction of neutral kaons was restricted to transverse K° momenta pt > 1 GeV, 
although the theoretical predictions should be more reliable for transverse momenta above 
3 GeV [0]. A comparison of data with this calculation is shown in Fig. [3]. In the high p t 
region the agreement is satisfactory. Below 2 GeV the calculation deviates significantly 
from the data. 



4.2 Comparison to Deep Inelastic Scattering and pp 

A comparison of K° production rates in jp and DIS is shown in Fig. |] and a similar 
comparison of A production rates in Fig. |5] for the non-diffractive event samples defined 
in section 3.2. The ^p and DIS rates are compatible in y* and at low pt- At high p t , 
where possible differences in the underlying hard subprocesses may manifest themselves 
(charm production via photon-gluon fusion, hard gluon radiation in DIS, resolved photon 
contribution in jp), the DIS data lie somewhat above the 7p data. 

It is of interest to investigate if such an agreement in rates extends to other processes, 
such as, for instance, pp collisions. Figure |6| shows a comparison of K° and A production 
in pp collisions at 200 GeV, obtained by the UA5 collaboration and the present photopro- 
duction measurement. The UA5 Kg data have been multiplied by 2 to convert the 
K s rates to K° and K° rates, and divided by the width of the rapidity range in which 
the measurement has been performed (\y*\ < 2.5). In the case of the UA5 A data the 
rapidity range is restricted to \y*\ < 2. The trigger of the UA5 experiment rejected events 
in which either only the p or p dissociate, as well as elastic events. This sample, termed 
non-single diffractive, may be compared to the non-diffractive samples of the present mea- 
surement as low-mass diffractive events are removed in both cases. It can be seen that the 
K° production rates (Fig. |6|a) are in good agreement. In contrast, the A rates in the photon 
fragmentation region in photoproduction are lower than in pp interactions. This may be 
attributable to the presence of an initial baryon in pp in this hemisphere from the incoming 
proton or antiproton. A difference between photon and target fragmentation region with 



respect to baryon production has also been observed in deep inelastic p,N scattering [TO 
Given the observed process independence of the low p t K° meson production rates, it 

may be asked if such a feature holds for other mesons. In order to answer this question, the 

rate of production of charged particles, which consist mostly of charged pions and kaons, 

was investigated in jp and DIS. 

For this analysis the charged particles are accepted in the extended rapidity range 0.5 < 

y* < 3.5. At the limits of this range the geometric acceptance decreases to 80%. 

Figure [7| shows the comparison of charged-particle rates observed in DIS and jp. Here 
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the data are extrapolated to pt = using samples of Monte Carlo events which include a 
full simulation of the effects of the HI detector. In Fig. [j]a, where all events contribute to 
the center of mass rapidity spectrum, it can be seen that the rate in DIS is significantly 
higher than in jp. When the (low mass) diffractive events are removed (Fig. [j]b), the rates 
are comparable in both samples within errors. 

The simple quark parton model does not predict an explicit Q 2 dependence of the 
particle production rates but it can be argued that perturbative QCD would introduce 
such a dependence PTfl . A small Q 2 dependence was reported by EMC |28| for charged 
particles, but this dependence decreases with increasing W. 

The present measurement does not indicate a significant change of non- diffractive pro- 
duction rates for K°, A and charged particles in the photon fragmentation region when Q 2 
changes from Q 2 ~ to Q 2 ~ 20 GeV 2 . This extends the observations made at higher Q 2 
by HI for charged particles and by ZEUS for strange particles [[29[ that there is little or 
no dependence of the average particle multiplicities on Q 2 . However, it should be kept in 
mind that the spectra investigated here are dominated by low p t particles. At higher p t , as 
mentioned earlier, the DIS and 7p kaon rates show a small discrepancy. This observation is 
consistent with results published by HI showing that for p t > 1 GeV the charged-particle 
rates increase with increasing Q 2 M. 



5 Conclusion 

The K° and A production cross sections have been measured in photoproduction and 
DIS with the HI detector at HERA, at an average 7p center of mass energy of about 
200 GeV, as a function of the transverse momentum and the rapidity, in a kinematic 
range corresponding to the photon fragmentation region. No significant enhancement of 
the non-diffractive K° and A production rates is observed when going from Q 2 ~ to 
Q 2 ~ 20 GeV 2 . This feature has been found to hold also for the charged-particle non- 
diffractive production rates. At higher p t , however, the K° rate in DIS is somewhat higher 
than in jp. The K° rate has been found to be similar to that measured in pp interactions 
at the same center of mass energy. However the measured A rate is lower than in pp. The 
measured K° spectra have been found to be broadly consistent with Pythia and NLO 
predictions, whereas the measured A spectra are significantly underestimated by Pythia. 
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> 2000 




7i + 7i"mass [GeV] 




1.1 1.11 1.12 1.13 1.14 



p7i mass [GeV] 

Figure 1: Invariant mass distributions for K s (a) and A (b) candidates in photoproduction, 
after applying all the cuts described in the text. The Kg candidates are restricted to the 
kinematic range \rj\ < 1.3 and 0.5 < p t < 5 GeV; the A candidates are restricted to the 
kinematic range \r)\ < 1.3 and 0.6 < p t < 5 GeV. The curves represent the fits as described 
in the text. 
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Figure 2: Measured K° (a,b) and A (c,d) total ep cross sections in photoproduction as a 
function of p t and pseudorapidity r] compared to PYTHIA in the kinematic range \r]\ < 1.3 
and 0.5 < p t < 5 GeV. The inner vertical error bars indicate the statistical error only and 
the outer error bars the statistical and total systematic error added in quadrature. The 
labels E665 and DELPHI correspond to Jetset parameters favored by these experiments; 
the label A s = 0.3 denotes standard Jetset settings. 
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Figure 3: Measured K° total ep cross section in photoproduction as a function of p t 
compared to NLO calculations in the kinematic range \r)\ < 1.3 and 0.5 < p t < 5 GeV. 
The inner vertical error bars indicate the statistical error only and the outer error bars the 
statistical and total systematic error added in quadrature. The NLO curves correspond to 
different values of the renormalization and factorization scales fi, namely pt, 2p t and pt/2, 
where pt is the transverse momentum of the K°. 
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Figure 4: K° rate as a function of p 4 (a) and center of mass rapidity y* (b) in 7p and DIS 
for non-diffractive events. The data are integrated over 1.3 < y* < 2.8 in (a), and over 
p t > 0.5 GeV in (b). The inner vertical error bars indicate the statistical error only and 
the outer error bars the statistical and total systematic error added in quadrature. 
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Figure 5: A rate as a function of center of mass rapidity y* in 7J9 and DIS for non-diffractive 
events. The data are integrated over p t > 0.6 GeV. The inner vertical error bars indicate 
the statistical error only and the outer error bars the statistical and total systematic error 
added in quadrature. 
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Figure 6: K° rate (a) and A rate (b) as a function of p t in 7p and pp (UA5) for non- 
diffractive events. The HI data are integrated over 1.3 < y* < 2.8 and the UA5 data 
are normalized to the same rapidity range. The inner vertical error bars indicate the 
statistical error only and the outer error bars the statistical and total systematic error 
added in quadrature. 
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Figure 7: The charged-particle center of mass rapidity distribution in DIS and 7p for all 
events (a) and non-diffractive events (b). The data have been extrapolated down to p t = 0. 
The inner vertical error bars indicate the statistical error only and the outer error bars the 
statistical and total systematic error added in quadrature. 
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